Electron-induced rippling in graphene 
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We show that the interaction between flexural phonons, when corrected by the exchange of 
electron-hole excitations, may place the graphene sheet very close to a quantum critical point char- 
acterized by the strong suppression of the bending rigidity of the membrane. Ripples arise then 
due to spontaneous symmetry breaking, following a mechanism similar to that responsible for the 
condensation of the Higgs field in relativistic field theories. In the presence of membrane tensions, 
ripple condensation may be reinforced or suppressed depending on the sign of the tension, follow- 
ing a zero-temperature buckling transition in which the order parameter is given essentially by the 
square of the gradient of the flexural phonon field. 



Introduction. — The discovery of graphene [D-Qi the 
two-dimensional (2D) metallic crystal made of a carbon 
monolayer, has spurred a flurry of research activity due 
to its unique electronic properties and potential applica- 
tions [4]. Its mere existence came somewhat as a sur- 
prise, since it was widely thought that thermal fluctua- 
tions should inevitably destabilize atomic-thin crystals. 
Graphene, however, is not only stable, but surprisingly 
robust and stiff [5[ , crucially due to the strong and direc- 
tional a bonds holding the carbon atoms together. From 
a theory point of view, these bonds provide the fixed con- 
nectivity that protects graphene's long-range crystalline 
order against thermal dislocations, and the in-plane shear 
modulus that allows it to remain flat below the crumpling 
transition temperature [fj, |7|. 

A remarkable and unexpected property observed in 
this flat phase, moreover, is its strong tendency to de- 
velop ripples, or long wavelength modulations of the out- 
of-plane displacements Q . X-ray diffraction experiments 
on suspended graphene revealed that this rippling is not 
simply a thermal excitation of flexural phonons around 
a flat average, but a true broken symmetry of the sys- 
tem that freezes into a corrugated average configuration 
[§] . Ripples are expected to have a significant impact on 
electronic transport in graphene Q. 

The origin of graphene's rippling has been heavily de- 
bated. In exfoliated graphene, ripples are correlated to 
some extent with the irregularities of the substrate on 
which the graphene sheet is deposited [l(|. But there 
is also evidence that they may arise in part as an effect 
intrinsic to the 2D membrane It has been proposed, 
for example, that variable length cr-bonds characteristic 
of carbon[ll|. or sp 3 orbital rehybridization driven by 
electron-phonon coupling [l2| could underlie ripple for- 
mation. 

In this work we argue that 7r-electrons in graphene play 
a crucial role in the ripple formation mechanism, since 
their coupling to flexural modes may lead to a strong 
renormalization of the elastic parameters. Above certain 
value of the electron-phonon coupling, we find that the 
system gets very close to a quantum critical point charac- 



terized by the strong suppression of the bending rigidity. 
We propose that, in this situation, ripples arise naturally 
by spontaneous symmetry breaking, following a mecha- 
nism similar to that responsible for the condensation of 
the Higgs field in relativistic field theories. We are able 
to characterize in this way a buckling transition that is 
present at the level of the zero-temperature field theory, 
with an order parameter controlled by the sign of the ten- 
sion in the membrane, and leading to a predicted buckled 
phase consistent with experimental observations. 

Softening of flexural phonons. — We describe the elas- 
tic deformations of a graphene sheet by the vector field 
u = (ui,U2,h), where ui,U2 represent the in-plane dis- 
placement with respect to the equilibrium position and h 
is the out-of-plane shift. The elastic energy of the mem- 
brane is expressed in terms of the strain tensor 

Uij (r) = i (diUj + djUi + dih djh) (1) 

Coordinate r runs over the 2D graphene membrane. The 
bare parameters involved in the continuum elasticity 
model for graphene arc the mass density p, the bending 
rigidity (k w 1 eV) and the in-plane shear (/jr;9 eV/A ) 
and bulk (/i+Aw 12.6 eV/A ) moduli. The action for 
the phonon fields reads 



S u = ^ / didMp(W-*(V^r-2M<-A(^r) (2) 



Respect to other conventional insulating membranes, a 
novel effect in graphene is that the electronic carriers cou- 
ple to the displacement of the sheet, in such a way that 
it gives rise to another source of interaction between the 
phonon fields IH 14 1 . We will consider the effect of elec- 
trons from the n bands of graphene, represented by the 
field ^(r), on the renormalization of elastic parameters. 
The strongest coupling between electrons and phonons 
comes from the on-site deformation potential (l5j . 

Sc-ph = / dtd 2 r * t (r)*(r) + g^d^h) (3) 
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FIG. 1. Lowest order corrections to the self-energy of flexural 
phonons (represented by a curly line) arising from (a) the four- 
phonon interaction Ko and (b) the exchange of electron-hole 
excitations (represented by the bubble with fermion lines). 
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FIG. 2. (a) Effective momentum dependence of k obtained 
from Eq. J7J for k = 0.8 eV, K = 20.5eV/A 2 , and values of 
the deformation potential g equal (from top to bottom) to 0, 
21 and 23 eV. (b) Scaling of k in the RG approach, for bare 
values of g equal (from top to bottom) to 0, 27 and 30 eV (for 
which the scaling of Ko and Ki = Ki/q c is also displayed in 
the inset), a is the C-C distance. 



with general momentum dependent couplings SWout(q)- 
Rotational invariance implies that at small momenta 
5i n and <7 ut must coincide. A microscopic analysis 
shows however that they can be very different at large 
momenta, with a general behavior <7i n /out(q) = 5(1 — 
c in/out|q| 2 + • ■ -)0- Detailed electron-phonon calcula- 
tions have shown indeed that <7i n has very small values 
< 1 eV at the K point of the Brillouin zone[l6j|. 

The energy of the electronic excitations is proportional 
to the Fermi velocity vf and (for a given momentum) 
much larger than that of other modes in the problem. 
Thus, we can first integrate out the electron degrees of 
freedom, which leads to a phonon interaction due to the 
exchange of electron-hole pairs, proportional to the static 
charge susceptibility %(q, 0) ' 



m\/Vp [17|. One may 
then integrate out the in-plane phonons, and arrive at 
the action for the flexural modes 



d 2 q dui 2 4 

K(q)u(q,u))u(-q,-uj) (4) 



d 2 q dui 
(2tt) 2 2^ 



FIG. 3. (a) Exchange of flexural phonons (curly lines) leading 
to logarithmic cutoff-dependence of the interactions Ko,Ki 
and K%. (b) Generic form of the diagrams building up the 
power series of the a field (dashed lines) in the effective po- 
tential, of which (a) is the first term. 



where u(q, ui) stands for the Fourier transform of 
^Pijdihdjh, Pij being the transverse projector [6], and 

K(q) = 2/i + A - . 9out |q|/, F - - - - _ (5) 

An electron-induced bulk instability occurs if the pole 
in the coupling function K(q) falls inside the first Bril- 
louin zone, leading to a vanishing velocity of in-plane 
longitudinal phonons. The possibility that such an in- 
stability could explain ripple formation [l3j| is ruled out 
however in a real graphene system, since experimentally 
the in-plane integrity of rippled graphene is intact. But 
even if the coupling <?i n (q) is not strong enough to desta- 
bilize the bulk, exchange of particle-hole pairs can make 
the flexural phonon coupling function K(q) negative (i.e. 
attractive), which indicates a potential instability of out- 
of-plane displacements. The effect of the negative terms 
in Eq. ([5]) is only significant at short wavelengths, and 
vanishes in the limit q — > 0. Then, in order to study the 
competition between positive and negative couplings in 
Eq. ([5]), we may approximate this coupling function by 
the constant term Kq and the dominant powers of |q| 



K(q) &K - Ki |q| - K 2 |qf 



(6) 



with K\ oc g 2 and Ki oc g 4 . 

The resulting interactions tend to enhance or suppress 
the low-energy rigidity k in opposite directions depending 
on their repulsive or attractive character. This is encoded 
into the corrections to the phonon self-energy S(q, ui) 
represented in Fig. [TJ Computing the full propagator 
of the flexural phonons from the expression D _1 (q, ui) — 

— ( K o//°)q 4 — S(q, ui), we find that the dressed bending 
rigidity «(q) is given by the self-consistent equation 



«(q) = Kq + 



1 

8^ 



d 2 p sin 4 ( 



ff(|q-p| 



q-pr 



VW) 
(7) 

This integral equation can be solved in different 
regimes of the couplings, see Fig. [2ja). We observe that, 
for a sufficiently large value g* of the deformation po- 
tential, /c(q) vanishes at a certain momentum, bouncing 
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back for smaller values of |q|. In this picture, g* plays the 
role of critical coupling above which we would find neg- 
ative values of the bending rigidity, suggesting a flexural 
instability in the system. In this situation, a renormal- 
ization group (RG) analysis may be applied in order to 
also capture the low-energy scaling of the couplings K n , 
which is quite significant when k becomes very small. 

The RG scheme proceeds by progressive integration of 
energy shells, starting from a high-energy cutoff E c = 
\J Ko/p Qc j ^0 approach the low-energy regime. The bare 
couplings K n are then corrected to lowest order by the 
exchange of two flexural phonons given by the diagram in 
Fig. EJa), which shows a logarithmic dependence on the 
cutoff q c . This can be absorbed into a renormalization 
of the effective couplings, leading to the scaling with the 
running cutoff q c — » 



dKn 
dq c 



n K„ 



3 
64^r 



(8) 



i+j 



A scaling equation can be also written for the rigid- 
ity k, as the self-energy diagrams in Fig. [1] display a 
logarithmic dependence on the cutoff q c , 



8k 
dq c 



3 K, 



16tt 



/ pn 



3 
16tt 



E 

n=l,2 



pn 



(9) 



The behavior of k obtained from the scaling equa- 
tions is shown in Fig. HJb) for different bare values 
of g. For small values of this coupling, we see that 
the effective rigidity has a steady increasing trend as 
q c — > 0. There is however a regime, corresponding to 
g ~ 20 — 30 eV, in which k is significantly reduced due to 
the electron-phonon interactions encoded into K 12 - For 
sufficiently strong g, k will be suppressed by several or- 
ders of magnitude, eventually driving the system into a 
strong-coupling regime K ^> \J pn 3 , and then bouncing 
back to become trapped in a range of very small values. 
As in the self-consistent approach, this points to a transi- 
tion into a new regime in which flexural modes are desta- 
bilized. Although the deformation potential is subject to 
a si gnifi cant uncertainty in graphene, with g ~ 4 — 30 eV 

mil, 

the bare values for the RG approach should cor- 
respond to the large estimates in that range, which apply 
to the interaction with high-energy electronic states [19(. 
We note that the possibility of such an electron-mediated 
rippling instability is consistent with the observed ab- 
sence of ripples in fluorinated graphene [23j, which differs 
with respect to clean graphene in the lack of conduction 
electrons from the 7r band. 

Spontaneous symmetry breaking. — Close to vanishing 
k, the renormalized action for h takes the form of a crit- 
ical scalar field theory with repulsive interaction Kq and 
a reduced cutoff E c . The ground state properties of this 
theory are nontrivial due to the role of quantum correc- 
tions, which may destabilize the classical (flat) solution 
leading to a ground state with broken symmetry. The 



standard technique to compute the equilibrium configu- 
ration of the s yste m is the minimization of the effective 
action S e s(h) [2J]. While symmetry breaking has been 
described before in the statistical mechanics of polymer- 
ized membranes [2 5], here we deal with the computation 
of the effective action of the quantum theory at zero tem- 
perature. This can be accomplished by introducing first 
an auxiliary field a to decouple the four-phonon interac- 
tion in (QJ, allowing to express the interaction term as 



S\nt — — 



2 J (2tt) 2 2tt 



cr(q,w) cr(-q, 



'Kou(-q,-u) 



(10) 
h into an 
Quantum 
what can 



The effective action is found by decomposing 
average field h av and quantum fluctuations h. 
corrections are obtained by integrating over h, 
be accomplished exactly in the formal limit of a large 
number of dimensions d of the ambient space containing 
the membrane. 

At large d, the effective action is built from diagrams 
with just one loop of fluctuating h fields, as represented 
in Fig. EJb). In the case of static field configurations 
M av (q,w) = <S(w)u av (q) and £r(q,w) = 5(w)cr(q), we get 
the contribution to the effective action 



iS. 



oo 

(1) _ \^ 



2 



1 K 



n/2 



n (27r) r 



|p|>A 



qd<A i=1 

d 2 p dujf 
(2tt) 2 2tt 



(2tt 



,-,27) 



sin 



>Jp,q 



i) 



(pu$ - np A ) 



~ (11) 



We focus here on the field configurations as A — > and 
the momentum of <r(q) goes to zero, since this is the 
regime where terms with higher powers of the field be- 
come increasingly infrared divergent. 

The perturbativc scries has to be summed first in (jlip 
to obtain a sensible result in the limit A — > 0. After 
factoring out the volume of the space, we obtain the con- 
tribution to the effective potential 



(i) _ 



err 



d 2 p dujp 
(2tt) 2 ~2V 

log (l 



2tt 



-0o- 



-pa; 2 . 



Kp 1 



2tt 



-0"o- 



Kp H 



(12) 



where cto = (A/27r) 2 er(q — > 0). Performing now the loop 
integral, we get from the sum of the zero-loop and the 
one-loop effective potential 



^eff 



1 

8^ 



-a + 2\J K o- 



1(2*) 



+ 



1 Kg 2 / / VA^Q CTQ 

8tt 2 16^ v /^ 3/2fT ° ^ ° S \^^pRE c 



2 )(13) 



Quite remarkably, Eq. (TTBj) exactly reproduces the 
structure of the effective potential for a relativistic scalar 
field theory in 3+1 dimensions 24j. This means that our 
model actually follows the same mechanism of symmetry 
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breaking characteristic of a Higgs field in particle physics. 
In close analogy with the analysis of Ref. [26j, the slight- 
est tension in our model will add to the potential a term 
proportional to 7M a v(q —> 0), destabilizing the minimum 
at the origin for 7 < 0. The mode cto can be integrated 
by the saddle-point method, showing that V c s will get 
then the typical "mexican-hat" shape as a function of 
w av , implying a minimum with V/i av 7^ 0. 

We find therefore that there is a buckling transition 
between the regimes of positive and negative tension in 
graphene. The nature of the critical (tensionless) theory 
is less clear, as there is still controversy about whether a 
scalar field theory may undergo symmetry breaking start- 
ing with a massless scalar field. Anyhow, in a typical ex- 
perimental setup, graphene is held in place by a scaffold 
or a substrate, which induces certain amount of tension 
in the membrane. This has been estimated to be of the 
order of ~ 1 eV / nm 2 [l2[ • In our analysis of the symme- 
try breaking, a negative tension 7 shifts the minimum of 
the effective potential to (A/27r) 2 u av (q — > 0) = — / -f/2K - 
This quantity turns out to be of the order of ~ 10 with 
the above estimate. Assuming that (A/2-7r) 2 scales in in- 
verse proportion to the volume of the space, we obtain 
that the average values of |Vh a v| must be of the order of 
~ 10~ 2 , which is consistent with the typical aspect ratio 
of ripples in graphene. 

In the tensionless limit, our model also shares with 
massless scalar field theories in 3+1 dimensions the in- 
triguing feature that the effective potential becomes com- 
plex beyond a certain value of the field. This kind of 
instability has led to speculate about the possibility that 
the absolute minimum of the effective potential may be 
away from the origin in the critical theory 26] . We hope 
that the connection between buckling and the Higgs con- 
densation may shed light into this question, while carry- 
ing out experiments at the low-energy scale of graphene. 

Conclusion. — We have seen that the interaction be- 
tween flexural phonons mediated by particle-hole excita- 
tions may place the graphene sheet very close to a quan- 
tum critical point characterized by the strong suppression 
of the bending rigidity of the membrane. The same ef- 
fect can be expected from other degrees of freedom which 
couple to the lattice deformations, such as absorbed im- 
purities. The effective potential of the zero-temperature 
theory displays then a mechanism of symmetry breaking 
that parallels that of the Higgs field in particle physics, 
the role of order parameter being played in graphene by 
the square of the gradient of the flexural phonon field. 
We have found that the system must be unstable to- 
wards a buckling transition that is the analogue of Higgs 
condensation, showing another way to employ graphene 
as a test ground of fundamental concepts in theoretical 
physics. 
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